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Abstract 

We detect an omission in the paper "Conditionally exactly soluble class of 
quantum potentials" by A. de Souza Dutra [Phys. Rev. A 47 (1993) R2435]. 
There, two strongly singular s— wave bound state problems have been claimed 
completely solvable in closed form. Unfortunately, all the displayed wave func- 
tions represented merely asymptotically correct (so called "Jost") solutions 
and did not satisfy the appropriate threshold boundary condition. We show 
that the incorporation of its standard form only leads to a very partial exact 
solvability at a single energy and for special couplings. 



For the two strongly singular s— wave potentials given, in the units h — 2/j, — 1, 
by the formulae 



and 



W = ^ + ^ + ^, G = 9 „=-A (2) 



A. de Souza Dutra [1] offered the explicit elementary wave functions as well as 
closed formulae for all their bound-state energies. One of the three couplings 
is not free: This entitled him to coin their "conditionally" exactly soluble 
(CES) status. In what follows we intend to demonstrate that in the sense of 
the Ushveridze's monograph [2] both these forces Vi^fV) only remain partially 
solvable at certain specific values of the energies E and couplings B. 

Our present main point is that all the solutions presented in ref. [1] still 
have to satisfy an appropriate and, for reasons to be made understandable 
here, forgotten boundary condition in the origin. Indeed, it is well known that 
for a central potential, Schrodinger equation — A\l/(f) + l / (|r|) , I'(r) = E^(r) 
degenerates to an infinite set of the ordinary (often called radial) decoupled 
differential equations 

- ^(r) + ^^W) + V(r)^(r) = E^(r), £ = 0, 1, . . . (3) 

for the separate angular-momentum components of the whole original wave 
function. The Newton's excellent review [3] summarizes the details. Under the 
assumption of the analyticity ofV(r) in the origin it shows that and why the 
standard physical requirement of normalizability of bound states ||\l/(r*)|| < oo 
is strictly equivalent to the integrability of their partial waves, 

ip{r ) € L 2 (0, oo). (4) 

For £ = 1,2, .. . the unphysical component ipirreguiarij') ~ r~~ e of the general 
threshold solution of eq. (3) is manifestly non-integrable near r 0. 
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In the s— wave with t = a more subtle argumentation is needed [3]. In 
practice, the subtlety is usually avoided by the replacement of eq. (4) by the 
boundary condition 



Even when we solve the ordinary harmonic oscillator the latter boundary condi- 
tion in the origin offers a more straightforward recipe for numerical calculations. 
Let us repeat: for analytic potentials, eqs. (4) and (5) are equivalent but the 
proof [3] of their equivalence immediately fails for the "very next" non-analytic 
V(r) Gr~ 2 , r ps 0, say in the Kratzer's solvable phenomenological model [4] 
with G ^ etc. One must re-analyze the whole quantization procedure anew, 
even for harmonic oscillator at G — > [5]. 

For all the similar singular forces with the finite limit in eq. (3), 



we have to re-define the dominant singularity + 1) + G = £(£ + 1). The new 



parameter £ — y [£ + |J + G — ^ enters then the modified threshold solutions 
ipreguiar( r ) ~ r ' C+l an d ^ 'irregular ij') ~ The irregular one is eliminated as 
manifestly violating the normalizability (4) at £ > 1/2. 

The latter bound means G > 3/4 in s— wave with £ = 0. Below such a 
strength of repulsion the Hamiltonian ceases to be self- adjoint. The conclusion 
is strongly counter-intuitive. Mathematically, the problem is serious. First 
spotted and analyzed by Case [6], it means that at G < 3/4, the textbook 
quantization of the Kratzer-like singular models is not unique at all. A more 
detailed discussion may be found in the literature (cf., e.g., [7] or [8]). In its 
light, physics community currently accepts a unique way of quantization which 
is, mathematically speaking, a mere regularization. It is often supported by 
the various sufficiently robust ad hoc arguments (cf., e.g., [5] on pp. 157 and 
167orref. [9]). 

For our present purposes, in the physical language of textbook [10], the 
correct recipe may be formulated as follows. 



\imib(r) = 0. 

r^o y ' 



(5) 



G = lim r 2 V(r) ^ 
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• In the domain of a weak repulsion we distinguish between the physical 
^reguiarif) ~ r c+1 and unphysical ipirreguiar(r) ~ r~ c . As long as both 
of them remain normalizable, we impose an extra, stronger boundary 
condition in the origin, 

lim^(r) = 0, Ge (0,3/4). (6) 

It coincides with (5) but its mathematical meaning of a convenient choice 
of the most plausible self-adjoint extension is different. 

• In the domain of weak attraction, both solutions ip re guiarij') ~ r ' C+1 an d 
^ 'irregular (j') ~ r ' C are compatible with eq. (6). In a sensible physical 
theory which distinguishes between the two, the replacement of eq. (6) 
by an even stronger artificial constraint is needed, 

lim^(r)/v^ = 0, Ge (-1/4,0). (7) 

• Below the lower bound G < —1/4 one cannot prevent the spectrum from 
collapse by any means. Particles would definitely fall in the origin. 

We may summarize: In practice, bound state solutions of the Schrodinger 
differential eq. (3) may be constructed in two ways, namely, 

• [RS] as the regular solutions ipreguiarij') constrained by the asymptotic 
normalizability condition 

^regular (R) =0, R — > OO; 

• [JS] from the so called Jost solutions ipj os t(r), always exhibiting the 
square-integrable asymptotic decrease by definition. 

The former regular-solution approach [RS] proves useful within the framework 
of the standard Taylor series method [11] and in non-numerical context [12]. 
Schrodinger equation (3) becomes converted into the exactly solvable two-term 
recurrences at q = (harmonic oscillator), into the three-term recurrences at 
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q — 1 (sextic forces) etc [13]. Rather unexpectedly, for all the positive integers 
q — 1, 2, . . ., a few bound states may still appear in an exact polynomial (i.e., 
terminating Taylor-series) form. An explicit construction of these exceptional 
elementary states is based on the solution of the Magyari's nonlinear algebraic 
equations [14]. They determine a few energy levels exactly and restrict also 
the free variability of the available couplings. 

In the latter context, potentials Vi ;2 (r) exhibit a certain incomplete dy- 
namical symmetry and play an exceptional role as quasi-exactly solvable in a 
certain narrower sense (cf. ref. [2] for more details). This would make the 
ambitious conclusions of ref. [1], if they were all true, even more important. 

Their analysis must be based on the alternative option [JS] which requires 
the threshold boundary condition (5), (6) or (7) [7]. This is a core of our 
present message. For the particular forces (1) and (2) such an approach has 
already thoroughly been tested numerically in ref. [9]. The Liouvillean [15] 
change of variables r — > x = r const and ip(r) — >■ x const x(x) has been employed 
there. As long as it leaves the form of the Schrodinger equation unchanged, 
it reduces all the bound-state problems with forces of the type (1) and (2) to 
their "canonical" equivalents with polynomial potentials 

V T (x) =ar- 2 + br 2 + cr A + ... + yr Aq + zr Aq+2 , a > -1/4. (8) 

On this basis, we may easily deduce the leading-order solutions (near the origin) 
also for our singular potentials Vi ;2 (r) of eqs. (1) and (2), 

^Pl,regular{^ ~ T ^ , ^Pl,irregular{^ ~ T ^ , 

1p2,regular{r) ~ r 5/6 , ^irregular (?") ~ r 1/6 . 

This is to be compared with the Dutra's wave functions: Say, for potential 
Vi(r) we may quote equation Nr. (9) from [1], 



ij{ D \r) = C r 1/4 exp 



1 o2 f.A/2 & \ 2 



--P 2 ir 1 ' 2 - 
T V 2E 



[3 (r 



1/2 _ B_" 
IE, 



(9) 
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Its energies E = — /9 4 /4 are parametrized by (3 = f3 n and numbered by an 
integer n = 0, 1, . . .. The formula also contains a certain normalization con- 
stant C = C n and Hermite polynomials H n {x). We immediately detect an 
inconsistency of the latter two equations near the origin. 

Similar observation is also made for ip2(j') from equation Nr. (13) in [1]: 
None of the Dutra's wave functions satisfies the physical boundary condition 
(7). An explanation of this obvious misunderstanding is in fact not too difficult: 
The solutions were merely constrained by the too weak (though, in practice, 
much more frequently encountered) and, hence, inapplicable threshold con- 
dition (5). We may summarize that the inconsequent use of the boundary 
conditions would lead to a physically absurd spectrum covering the whole real 
line, E G (— oo, oo). 

The Dutra's "non-anonymous" (i.e., Hermite-polynomial) solutions have al- 
ready evoked a non-negligible response in the current literature. As an example 
one might quote the paper [16]. Its authors relied on the physical correctness of 
the Dutra's argumentation and were misguided in their mathematical appreci- 
ation of the role of supersymmetry in the similar problems. Still, the majority 
of their argument remains valid. Hence, let us show, in the conclusion, how 
the Dutra's exceptional solutions could be "saved" for similar applications. 

Obviously, one has to incorporate simply the necessary constraint (7). An 
inspection, say, of our sample eq. (9) reveals that ip[ D \r) satisfies condition (7) 
if and only if its Hermite-polynomial component acquires an exact nodal zero 
in the origin. In terms of the known numbers X — X(n, k) (calculated as the 
k— th nontrivial zeros of H n (X), cf. Table 1) this requirement, unfortunately, 
fixes the non-Coulombic coupling as a function of the energy E = — /9 4 /4, 

B = ±X0 3 ^O. (10) 

This makes both these values coupled to the additional (in fact, Magyari's [14]) 
constraint. As a cubic equation for the energy E it appears under Nr. (8) in 
ref. [1]. This algebraic selfconsistency condition must be combined with eq. 
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(10). The resulting polynomial equation in (3 (of twelfth degree!) is easily 
factorized in closed form. Its real roots we need are 



They all exist for any A < 0. This is an important conclusion: let us note that 
eq. Nr. (8) of ref. [1] re-appears as equation Nr. (16) in ref. [16], etc. 

For illustration, let us finally fix the scale A = — 1 and display the first 
few non- numerical specifications of energies E = —f3 A and their couplings (10) 
in Table 2. The same parameters are to be used also in the definition (9) of 
the correct bound-state wave function. Mutatis mutandis, the entirely parallel 
"return to validity" applies also to ip2 D \ r ) m [!]• We omit the details here, 
re-emphasizing only that both the Dutra's expressions V'iT( r ) are elementary 
and still satisfy the Schrodinger differential equation, exhibiting also the correct 
asymptotic behaviour. Thus, we may return, say, to the paper by Dutt et al 
[17], originally motivated by ref. [1] as well. In the light of our present notes, 
the importance of the latter paper increases: Its authors have, involuntarily, 
found and constructed the first CES example in one dimension! 
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Table 1. Non- vanishing zeros X — X(n, k) 
of the first few Hermite polynomials H n (X). 



k 


4 


3 


2 


1 


n 





















1 
2 




-VV2 


>/V2 




3 




-^3/2 


^3/2 




4 


-^(3 + VE)/2 


-^(3 - >/6)/2 




>/(3 + V6)/2 


5 


-y/(5 + VT0)/2 


-y/(5-ym)/2 


v/(5 - Vl0)/2 


y/(5 + y/W)/2 



Table 2. Parameters of the first few simplest quasi-exact states in Vi(r). 
M counts nodes in ipi(r): M = means ground state, etc. 



M n k 


fixed coupling B' — B/8 


binding energy E 


2 1 

3 1 

4 1 

5 1 


^/l/9 3 - 0.0370 
y^/ll 3 - 0.0475 
^(3 + y/6)/(15 - V6 f - 0.0525 
^(5 + v / 10)/(17- v^) 3 - 0.0555 


_(4/9) 2 - -0.197 

-(4/11) 2 0.055 

-[4/(15 + V6)] 2 0.029 

-[4/(17 + VlO)] 2 0.018 


1 2 2 

4 2 

5 2 


-y/l/9 3 - -0.037 
y/(3- a/6)/(15 + >/6) 3 -0.0102 
^(5- v / 10)/(17 + v/lO) 3 -0.0150 


_(4/9) 2 - -0.132 
-[4/(15 - VQ)} 2 - -0.047 
-[4/(17- v^)] 2 0.028 


2 3 2 
4 3 


-y/3/U 3 - -0.047 
- ^(3 - >/6)/(15 + v^) 3 - -0.010 


-(4/11) 2 0.055 

-[4/(15 - v^)] 2 - -0.047 


3 4 4 
5 3 


-^/(3 + >/6)/(15- v^) 3 - -0.053 
-y/(5 - v / 10)/(17+ v^) 3 0.015 


-[4/(15 + v^)] 2 - -0.029 
-[4/(17- v^)] 2 - -0.028 


4 5 4 


-^(5 + v / 10)/(17 - v^lO) 3 0.056 


-[4/(17 + VlO)] 2 0.039 
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